The Diversity Potential of Relay Selection with 
Practical Channel Estimation 

Diomidis S. Michalopoulos, Member, IEEE, Nestor D. Chatzidiamantis, Student Member, IEEE, 
Robert Schober, Fellow, IEEE, and George K. Karagiannidis, Senior Member, IEEE 

Abstract 

We investigate the diversity order of decode-and-forward relay selection in Nakagami-m fading, in cases 
where practical channel estimation techniques are applied. In this respect, we introduce a unified model for 
the imperfect channel estimates, where the effects of noise, time-varying channels, and feedback delays are 
jointly considered. Based on this model, the correlation between the actual and the estimated channel values, 
p, is expressed as a function of the signal-to-noise ratio (SNR), yielding closed-form expressions for the 
overall outage probability as a function of p. The resulting diversity order and power gain reveal a high 
dependence of the performance of relay selection on the high SNR behavior of p, thus shedding light onto 
the effect of channel estimation on the overall performance. It is shown that when the channel estimates 
are not frequently updated in applications involving time-varying channels, or when the amount of power 
allocated for channel estimation is not sufficiently high, the diversity potential of relay selection is severely 
degraded. 

In short, the main contribution of this paper lies in answering the following question: How fast should 
p tend to one, as the SNR tends to infinity, so that relay selection does not experience any diversity loss? 
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I. Introduction 

Wireless relaying technology has been recently proposed as a method that promises significant performance 
improvement in wireless communications without any power increase Q, 0. Among the most common 
relaying techniques is so-called relay selection which has been extensively analyzed in the literature |[3l-f7l. 
In relay selection, the system is able to select a single relay out of the set of available relays, in order 
to take advantage of the multiple paths available and thus achieve spatial diversity. It has been shown that 
activating only the relay with the strongest instantaneous end-to-end channel represents a bandwidth-efficient 
alternative to all-participate relaying, since, on the one hand, the same diversity order is achieved, yet on 
the other hand, the excessive bandwidth usage that the activation of multiple relays entails is avoided Q. 

Most of the literature dealing with relay selection in fading channels has assumed that perfect channel 
state information (CSI) is available at the terminal where the decision on which relay to activate is made. 
This assumption, however, may not be true in practical scenarios where the channel changes rapidly enough, 
so that the CSI available at the selecting terminal is outdated. In addition, if the power allocated to the pilot 
symbols is not sufficiently high, the noisy channel estimates may lead to suboptimal relay selection. The 
above two cases reveal the vulnerability of relay selection to imperfect channel estimation, and constitute 
the main rationale for conducting a thorough outage and diversity analysis of relay selection in scenarios 
with imperfect CSI in this work. 

In fact, the case of relay selection under outdated CSI and Rayleigh fading has been recently studied in [ 8 ]- 
IPT5TI . where interesting results on the outage probability and diversity order were derived. Nonetheless, these 
works consider only a special case, since they assume that the CSI imperfection stems only from delayed 
feedback. This may not always be the case in practice, since channel estimates may also be impaired by 
time- varying fading and channel noise. In the very recent works lfT6ll , IfTTT . the effect of noisy channel 
estimates is also included in the performance analysis. However, these works are based on the assumption 
that the same estimates are used for both relay selection and detection, leading to zero diversity order. Such 
assumption may not always be true in practical scenarios where the number of pilots used for relay selection 
and that used for symbol detection may not be equal to each other. 

In light of the above, the contributions of this paper are summarized as follows. 

• We conduct an outage analysis of relay selection with imperfect CSI, which is general enough so as to 
account for both the effects of noisy and outdated channel estimates, integrated into a unified model. In 
particular, a closed-form expression for the outage probability is derived, which incorporates all effects 
that can cause imperfect channel estimates in practical applications. The considered channel estimation 
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techniques include the cases of estimation in noisy static channels; estimation in noiseless time-varying 
channels, and estimation in noisy time-varying channel with the aid of finite impulse response (FIR) 
and infinite impulse response (IIR) channel prediction. 

• The amount of CSI imperfection is reflected by the correlation coefficient between the actual and the 
estimated channel values, p, which is modeled as a non-decreasing function of the signal-to-noise ratio 
(SNR). As a result, the asymptotic outage behavior of relay selection is determined by the speed of 
convergence of the correlation coefficient to unity, as the SNR approaches infinity. The diversity order 
of relay selection with imperfect CSI is thereby derived, shedding light onto the diversity loss caused 
by imperfect CSI, with implications for the design of channel estimation techniques. 

• In contrast to other relevant works in the literature, where Rayleigh fading channels were assumed, 
the versatile scenario of Nakagami-m fading is studied. It is shown that the resulting diversity order is 
directly proportional to the fading shape parameter, m. 

Overall, the general conclusion of this paper is that the level of CSI imperfection plays an important 
role in the overall performance of relay selection, considerably affecting its diversity potential. A detailed 
discussion on the diversity order of relay selection for several practical channel estimation techniques is 
presented in Section IVI-Q and corresponding numerical examples are given in Section IVIII These results 
are based on the exact outage analysis conducted in Section [TV] for Rayleigh fading and certain channel 
estimation techniques, and extended to Nakagami-m fading in Section [V] Prior to the outage analysis, the 
unified model that incorporates the effects of noisy and outdated channel estimates is presented in Section 
iriTl The system model of decode-and-forward (DF) relay selection with imperfect CSI is given next, in 
Section |nj 

II. System Model 

Let us consider a cooperative relaying system which consists of a single source terminal, S, N DF relays 
which are denoted by R{, i = 1,...,N, and operate in the half-duplex mode [U, and a single destination 
terminal, D. 

Channel Model: Let Iiab denote the complex channel between nodes A and B, where A,B£ {S, D, R4 : 
i = 1, N}. Moreover, Rayleigh distributed fading in each of the participating links is assumed, implying 
that Iiab is complex Gaussian random variable (RV). The versatile scenario of Nakagami-m distribution of 
fading is considered in Section [V] In addition, since in this work we focus our attention on the asymptotic 
properties of relay selection under imperfect CSI, we assume independent and identically distributed (i.i.d.) 
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fading in each of the links involved. Moreover, the fading is considered slow enough such that Kab remains 
constant during the transmission of one frame. 

Let jab represent the instantaneous SNR of the link between terminals A and B, i.e., jab = \^ab\ 2 /Nq, 
where A^o is the additive white Gaussian noise (AWGN) power. Due to the i.i.d. fading assumption, the 
average SNR in each of the links involved is identical, and denoted by 7. Moreover, we use the notation 
fx(-) and Fx (•) to refer to the probability density function (pdf) and the cumulative distribution function 
(cdf) of RV X, respectively. 

Relay Selection Process: Among the available relays, only a single relay is activated in each transmission 
session, based on the selection cooperation protocol (H. In particular, the relay selection procedure is 
completed in two phases, as follows. In the first phase, the relays that can successfully decode the message 
form the so-called decoding set, denoted by S. 

Mathematically speaking, the decoding set S is defined as S = {R{ : 75^ > T}, where T denotes the 
outage threshold SNR, defined as the maximum SNR value that allows decoding; T is related to the target 
data rate, r, through T = 2 2r — 1. In the second phase, the destination collects the estimated CSI of the 
Ri-D links with i : Ri 6 S, and activates the relay with the strongest Ri-D channel. 

A fundamental principle throughout this paper is the fact that the relay selection is based not on the actual 
channel values but on their estimates, which are generally not equal to each other. In this respect, let Kab 
denote the estimate of channel Iiab, so that jab represents the estimated value of jab, as seen by the 
destination. Hence, denoting the selected relay by R K , we have 

K = arg max j^d- (1) 

The CSI imperfection is assumed to affect the relay selection process, but not the symbol detection at the 
destination. This is because the number of pilot symbols used for detection is typically higher than that 
used for relay selection, and the channel estimates for detection can be updated more frequently. A detailed 
description of the considered imperfect CSI model follows. 

III. Imperfect CSI Model 

The physical causes of the considered CSI degradation are the time-varying nature of fading channels, as 
well as finite pilot symbol power. In this work, both of these causes of imperfect CSI are integrated into a 
unified model, as shown below. 

The level of CSI imperfection is quantified by the correlation coefficient between the actual squared channel 
envelope value, |/iab| 2 , and its corresponding estimate, |/iab| 2 , where A and B can be any terminals of the 
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set {S,D,Ri : i = 1, ...,iV}. This coefficient is denned as (in the sequel, all channel indices are dropped 
due to the i.i.d. fading assumption) 

_ e ((\h\ 2 - n h ) (H 2 -^)) 
P a \h\ 2(J \h\i 

where £lh = E^\h\ 2 ^, Q~ h = E (\h\ 2 ^j with E (■) denoting expectation, and ax denotes the standard 
deviation of RV X. 

Note that p reflects the effect of imperfect CSI on the SNR of the selected relay and hence on the overall 
performance of relay selection with imperfect CSI. For this reason, the subsequent analysis focuses on 
expressing the performance degradation as a function of p, so that all physical phenomena that cause CSI 
degradation are, in fact, incorporated into p. 

A. Versatile Imperfect CSI Case 

In accordance with the intuition that channel estimation is performed in noisy environments, let us consider 
the versatile scenario where p is a function of the SNR. That is, 

1 - P = 9 (7) (3) 

where g (•) is generally a non-increasing function of its argument, with < g (7) < 1. In order to obtain 
insight into the asymptotic dependence of the CSI error on the SNR, we expand g (7) into a Puiseux series 
llllfl . so that for high SNR we have 

g (7) = bT a + o {T a ) (4) 

where a and b are positive constants with < b < 1, and o (7"°) is defined such that Hindoo o (7~ a ) / (7 _a ) = 
0. 

It is emphasized that the versatile imperfect CSI model considered in © and (0]) is general enough 
to accommodate the cases of imperfect CSI due to noise impairment and the time-varying nature of the 
underlying channels. Next, we study the scenarios of CSI imperfection in static and time-varying Rayleigh 
fading channels, separately. 

B. Static Channels, Noisy CSI 

Let us assume the case of static channels, where channel estimation is implemented via averaging over L 
noisy pilot symbols. As a result, g (7) in (f3]) is a decreasing function of 7, for which lim^oo g (7) = holds. 
Moreover, let us consider the scenario where the power allocated to pilot symbols, £ p , is not necessarily 
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equal to the power allocated to data transmission, £d- We allow the ratio of £ p over to be SNR-dependent, 
so that 

£ p = PT£d (5) 

where (3 is a positive constant and a is a constant, the sign of which determines whether £ p increases 
or decreases with SNR. The estimated channel values are expressed as h = h + n p , where h and n p 
denote the true channel component and the remaining noise component, respectively. Given that the channel 
estimates are derived by averaging over L pilot symbols, the noise variance of the estimation process equals 

al p =N /(£ p L). 

Lemma 1: The correlation coefficient, p, between \h\ 2 and \h\ 2 , is given by 

sh _ n h _ L?r +1 _ 
9 n h n h + ai + i {) 

Proof: Since h equals a linear combination of complex Gaussian RVs, \h\ 2 is exponentially distributed. 
Hence, it follows from the theory of the moments of exponential RVs that E (\h\ A ^ = 20^; E (|N 4 ^) = 2fi|. 
Using this result, the proof follows from (O after algebraic manipulations, in conjunction with (f5]) and the 
fact that 7 = £ d ft h /N . ■ 
Expanding ([6]) in a Taylor series for 7 — > 00 and using (O, we obtain for a > — iQ 

9 (7) = 1 - p = jzT {a+1) + o (r (a+1) ) ■ (v) 

Therefore, using (@]), from (O we have a = a + 1; 6 = 1/ (PL). 



C. Time-Varying Channels 

Next, the case of time-varying Rayleigh fading is studied, where the maximum Doppler frequency on each 
of the participating links is assumed identical, and denoted by Moreover, the autocorrelation function 
of the complex channel h is denoted by ph{Td); based on the Jakes' model |19], ph{Td) is given as 
Ph {Td) = ^hJoi^fdTd), where Jo (•) denotes the zeroth order Bessel function of the first kind [20, Eq. 
(8.411)]. 

1 ) FIR Channel Prediction: In time- varying environments, the channel estimation can be improved by 
utilizing the CSI available from previous time instances, so that the channel estimates are derived through 
a channel prediction process l2Tl . Let us consider an FIR channel prediction filter of length L, and denote 
the time interval between consecutive CSI acquisitions by Td. In such case, following the analysis in |2T1 . 

'The case of a < — 1 yields p = for 7 — > 00, and is out of the scope of this paper. 
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the predictor coefficients can be optimized so as to yield the minimum squared error between the actual and 
the predicted channel values, a\, resulting in 

a 2 e =n h -^R~\. (8) 

In ([8]), Uh denotes the L-dimensional autocorrelation vector, i.e., u^ = [ph (—To) , Ph {—LTd)] T , R denotes 

an L x L symmetric Toeplitz matrix, the first row of which is given by [ph (0) + Nq/£ p , ph (—Id) , 

Ph (—LTd + Td)], and (•) denotes the Hermitian operator. Hence, p is derived by combining (0 and ([8]), 

as 

p = u»R- l u h /n h . (9) 

It follows from © that the CSI error can be expressed as a function of the SNR as g (7) = 1 — p = 
l—u^R~ 1 Uh/^lh- Interestingly, it is noted that in the high-SNR regime and for fd>0,g (7) converges to 
a finite non-zero constant, i.e., 

Jim 0(7) = 1 -uf R- 1 ] u h /fl h = b>0 (10) 

implying that the CSI error is independent of the SNR. Hence, considering (O, it follows that for the case 
where the channel estimates are obtained through FIR channel prediction, a = holds. 

Ideal but Outdated CSI: This special case of channel estimation was considered in Il8l- lfl4l . and in 
fact corresponds to noiseless FIR channel prediction with a one-tap predictor (L = 1), and is dubbed as 
"outdated CSI" here. It implies that the CSI based on which the "best" relay is selected is noise-free, yet 
the selection of the "best" relay is not based on the current time instant but on a previous one, because of, 
e.g., a feedback delay. Based on (©, it can be shown that the correlation coefficient, p, for the outdated CSI 
case equals p = p\ (— T^) /Q"^, a result which is in accordance with (SI, J9J. Moreover, g (7) is a constant 
function in this case, and thus (0]) yields a = 0; b = 1 — p\ {—Td) /O?. 

2) IIR Channel Prediction: Let us now extend the channel prediction case to the scenario where the 
number of pilot symbols participating in the prediction process are infinitely large. As shown in Appendix 
lAl this scenario leads to a correlation coefficient of 

P = 1 - exp (r d j f _ A In [S hh (e^) + (/37 Q+1 )^] 4f) (/Jf^T^^ + ^f^T' dD 

where Shh (•) represents the Fourier transform of ph (■). Hence, combining ([3]) and (fTTT) . we obtain for high 
SNR 

u - (a + l)(l-2f d T d \ 

g (7) = exp (r d J ^ In [s hh (f 2 ^)] dfj /H 1 " 2 ^^ ' : '• (12) 
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Consequently, it is concluded that the parameters a and b of the asymptotic dependence of the CSI error on 
the SNR are given by a = (a + 1) (1 - 2f d T d ) and b = exp 0r d f[ d fd In [S hh (e j2irfTd )] dA p~^~V^ d )_ 

The reader is referred to Table U for an overview of how the parameters a and b are derived for the practical 
channel estimation scenarios considered in this paper. An asymptotic performance analysis of suboptimal 
relay selection follows. 

IV. Outage Analysis of Relay Selection with Imperfect CSI in Rayleigh Fading 

The outage probability is defined as the probability that the overall SNR lies below a given threshold, 
denoted here by T, i.e., P out = Pr{7 K < T}, where k denotes the index of the selected relay and j K is the 
corresponding end-to-end SNR. Observing that for all channel estimation scenarios considered in Section 
Hm h is obtained as a linear combination of complex Gaussian RVs, it follows that h is also a complex 
Gaussian RV. Hence, 7 is exponentially distributed. Consequently, the conditional pdf of the actual SNR, 7, 
conditioned on its estimate, 7, is obtained from [22, Eq. (2.11)] as 



/ 7 |7 ( x \y) = — T o 2 " — : r-= (13) 

where 7 denotes the average estimated SNR and Iq (•) denotes the zeroth order modified Bessel function of 
the first kind |20, Eq. (8.447.1)]. It is emphasized that since the parameters 7 and 7 are not necessarily equal 
to each other, which is in contrast to the outdated CSI case treated in flU, O, the diversity investigation 
of relay selection with imperfect CSI under the general imperfect CSI assumption requires that we conduct 
a new outage analysis for our scheme. This outage analysis is similar to that in |8], yet the corresponding 
expression for the (x\y) used in flHJ is substituted by (fT3l) . 

In particular, based upon the mode of operation of the selection cooperation the outage probability is 
expressed as 

N 

P out = Pr{S = $} + Y J F lRKD (T\ \S\=l)Pr{\S\=l} (14) 
1=1 

where |<S| denotes the cardinality of S and is the empty set. Because of the i.i.d. assumption for the fading 
in the S-Ri links, the second term within the sum in (fT4l) is given by [8, Eq. (7)] 

/ 



Pr{|5|=/} 



T 

1 — exp I — — 



v -' ' IT\ 



cxp(-— y (is) 



7, 

Furthermore, by defining Ai as the event that the ith relay out of / relays is selected, i.e., A4 : i = k, 
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F 1RkD (T I |<S| = I) is expressed as 

l 

F 1RkD (T\\S\=1) = ^F 7BreD (T| \S\=l, ^)Pr{A} 

i=i 

oo 

( S ll S 2)/7 RlD |A, (a^l^da^i (16) 

where we used the fact that Pr{Aj} = 1/Z, because of symmetry. The conditional density of 
conditioned on Ai is derived as 

/Vl^ (X2 ^ = Pr^I"} = Pr {^.| 11 F ^ fe) = 0^)^ (x 2 ) • (17) 

Consequently, substituting (TT3T > and ([T7l ) in (fT6l ) we obtain an expression for F lR D (T\ \S\ = I) which 
coincides with (8j Eq. (8)], where the case of outdated CSI was considered. This leads to an interesting 
observation which is summarized below. 

Under i.i.d. Rayleigh fading and assuming correlation coefficient p between the actual and the estimated 
SNR in each intermediate link, the outage probability of relay selection with imperfect CSI, P ou t, expressed 
as a function of p, is given by the same formula, irrespective of the channel estimation technique used. 
Equivalently, P ou t is independent of 7, a fact which can be explained by noting that only the relative values 
of % are relevant for relay selection, not their absolute values. Hence, scaling all the estimated SNRs by the 
same factor does not affect the relay selection process. Therefore, the outage probability of suboptimal DF 
relay selection for any channel estimation technique is as shown in [8 , Eq. (2)]. It is emphasized, however, that 
different channel estimation techniques lead to different dependences of p on the SNR, resulting ultimately 
in different diversity behaviors. The diversity order of relay selection with imperfect CSI will be studied in 
detail in Section [VI] 

V. Outage Analysis in Nakagami-tti Fading 

Let us now consider the case where the fading in all channels follows the Nakagami-m distribution |[23l . 
In this case, since the distribution of h is unknown for the unified imperfect CSI model, we confine ourselves 
to investigating the performance of relay selection with imperfect CSI for the three special cases presented 
below. 

A. Time-Varying Channels: Outdated CSI 

Recall from Section IIII-C 1 1 that this case corresponds to noiseless FIR channel prediction with one tap. 
Consequently, h represents a delayed version of h, and h follows the same distribution as h, so that the 
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joint pdf of 7 and 7 is obtained from the bivariate Gamma distribution ll24l . which is simplified using |[20l 
Eq. (9.210/1)] to 



00 (m);m 2m+2 i 



h, *a) = g ., { [ _ p)m+ 2 3 II f 



^ expj p^j- 



(m + j) 7 m +J 



(18) 



j=0 ^ fJ i=l 
where (x)^ denotes the Pochhamer symbol defined in |20l pp. xliii]. 

Following the same steps as in (fT4b- ([T7T ). and using the fact that, the pdf and cdf of the SNR for the 

Nakagami-m fading model are given by / 7 (7) = m m 7 m_1 / [7 m r(m)] exp (—7717/7) and F 7 (7) = 1 — 

T (m, ™y) /r(m), respectively, in conjunction with the binomial expansion ll20l Eq. (1.110)], we obtain 



the equivalent expression for (fT6l ). pertaining to Nakagami-m fading, as 



(TIP! = «) = «/ E-T7I 



2 1 00 



2m+2j 



m+ j — 1 
1 



^ -/o) m+2j 'r(m)r(m + j)7 2 ( m +J) 6XP V (1 - p)l 



l-i 

E 

i=0 



I - 1 



r (m) 



m+j-1 



exp 



mx2 



t(i-p) 



r ( m, —X2j dx2 j dx\. 

(19) 



It is observed from (fT9l that in order to derive a closed-form expression for F jR D (T | \S\ = I), the following 
integral needs to be solved 

roo 

I(fi,a,m,fi,j)= / x^ exp (— ax) V 3 (m, f3x) dx. (20) 
Jo 

For the case of m 6 Z, the integral in (l20b is evaluated as illustrated in Appendix iBl Hence, F 7s D (T | |5| = / ) 
can be derived, using (ETJ), (H Eq. (3.351/1)], and [20, Eq. (8.352/6)] as 

'Z - r 



_ ^ ^(r(m + j)-r(m + j, jzl 

F^ D (T\ \S\-l)-l^ jKl - p yT(m)T(m + j) ^ 



(-l)>(m + j>,i) (21) 



where we have set 



( 



r(* + i) 



1 



E 



50+?l+---+?m-l=« 



771—1 



v 



m— 1 



1 + E 



(22) 



9=0 



/ 



Consequently, a closed-form expression for the outage probability is obtained by combining (fl4l ). (fT5l ). and 
<l2Tb . yielding 



out 



r (m) 

^ (-ir 



00 JV l-l 

+ EEE' 

j=0 i=l i=0 



— jo)^T (m) r (m + j) 



r (m + j) — r I m + j, 



F[m, fT 
r (m) 

mT 



N-l 



(1 " P)7 



Z - 1 



r (m) 
^ (m + j, p, i) . (23) 
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It is noted that, for practical SNR values, i.e., 7 < 30dB, the infinite series in (1231) converges after a finite 
number of terms, not greater than 100. 

B. Time-Varying Channels: FIR Channel Prediction with Large L and IIR Channel Prediction 

In this case, the channel estimate h is obtained as the weighted sum of a large number of observations 
in time-varying fading scenarios (/^ > 0) ETTl . Hence, it follows from the central limit theorem that h is a 
complex Gaussian RV, so that \h\ 2 is exponentially distributed with average value denoted by 7. 

The joint pdf of 7 and 7 is obtained from ll24l Eq. (10)] by setting mi = m and 77J2 = 1, leading to 

exp U + — t f! m-l;ffl + i; _; \ ) (24) 



r (j + 1) r (j + m) \ 1 - p V 7 i J J \ ' '7(1-/5) 

where 1F1 (•;•;■) i s the confluent Hypergeometric function defined in [20, Eq. (9.210/1)]. Following the 
same procedure as in (fT4b- (fT7l) . we obtain the conditional cdf of jr k d as 



f, RkD (t\ \s\ = n = 



7 7 U-p, 



r _ <ti=0 [l + z(l-p)]' +1 r(m + j 

T 



exp(-^— 1 lP 1 ^-l; ro + J ;__jcix 1 (25) 

where p is given in (O and (fTTT) . It is observed that, for the same reasons addressed in Section ITVl 
F jR ^ D (T I |<S| = Z) is independent of 7. 

In case of m = 1, using (201 Eq. (3.351/1)], HB Eq. (8.352/2)] and the fact that 1 F 1 (0; 6; z) = 1, <|25]> 
reduces after some algebraic manipulations to 



i=0 v ' /J j=0 



(l-p) 7 l+i(l-p) 



t=0 t! eo-'W 



(26) 

As a cross check, it follows from [20 , Eq. (0.231)] and the infinite series representation of the exponential 
function |20, Eq. (1.211/1)], that (l26l) is equivalent to (U Eq. (8)]. In case of m > 1, (1251 ) in conjunction 
with (201 Eq. (3.351/1)] yields 

f <t\ 191 = /) = /vvv (YH-ira-pr^ r(m-i+jo 

111 j ^&fr [i + ^(i-p)F +1 r(m-i)r(m + J + /,)fc! 

x + j + fc) -T (m + j + k, ■ (27) 

The overall outage probability follows then from (1271 ) (or d26l ). if m = 1) and (fl4l . 
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C. Static Channels: Noisy CSI in the high SNR Regime 

Under the high SNR assumption, it is valid to assume that h = h + n p is also Nakagami-m distributed, 
with 7 = 7. Consequently, the outage probability for this scenario is given by (1231) . where p is given in the 
ensuing Lemma. 

Lemma 2: The correlation coefficient between \h\ and \h\ for Nakagami-m fading is given by 



L . (28) 



Proof: The proof is similar to that for Lemma [Qby using the fourth moment of a Nakagami-m distributed 



RV, E 



\h\ L 



m ft 



It is noted that (1281 ) reduces to © for m = 1. Expanding (1281 ) in a Taylor series for 7 — >• 00 and using 
©, we obtain for a > —I 

9 (7) = 1 - P « ^7" {Q+1) + o (7~ (a+1) ) • (29) 

Next, we shed light onto the asymptotic behavior of relay selection with imperfect CSI in Nakagami-m 
fading. 

VI. Diversity Analysis 

A. High-SNR Analysis 

Here, we present a high SNR outage expression, which is used as a stepping stone for deriving the 
diversity order of relay selection with imperfect CSI in Nakagami-m fading. For simplicity of exposition, 
this expression pertains to the cases of outdated CSI and noisy CSI for Nakagami-m fading, as well as 
FIR and IIR channel prediction in Rayleigh fading; an expression for FIR and IIR channel prediction in 
Nakagami-m fading follows, likewise, from d27T ). 

For sufficiently high values of 7, we have 



T { m >fy) = r(m)-7(m,f2/j ^ r (m) - = ^ _ ^m-y 

r (m) r (m) r (m) r (m) 7 m 

which holds based on the series representation of the lower incomplete Gamma function, / y(s,x) [20, Eq. 

(8.354/1)], in conjunction with the fact that as x — > 0, 7 (s, x) w x s /s. Therefore, by substituting (l30l in d23l) 

and setting p = 1 — &7 _a , as implied by (@]), we obtain an alternative expression for the outage probability 
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as a function of the parameters a and b for high SNR as follows 

( m m - l y m \ N 1(1) fl-bj- a \ n fm™- 1 ^ V_ ' 



T(m)j m J ^—^^—^j\r(m)T(m + j) \ a J \T(m)j r ' 
x(!-T>^J (r(m + j)-r^ + j ) ^JjG(m + i I a,6 > Z) (31) 



where we have set 

2-1 



G (m + j, a, 6, = E f • ^ ( m + 1 " 6 ^" a ' • 

8=0 \ Z / 



(32) 



B. Diversity Order 

An important result derived from the high SNR analysis of Section IVI-AI is the diversity gain of the 
scheme under consideration, which is summarized in the ensuing theorem. 

Theorem 1: The diversity gain of relay selection with imperfect CSI with practical channel estimation 
and Nakagami-m fading is given by 

m [a (N - 1) + 1] , if a < 1 



G d = < 

Proof: The proof is given in Appendix O 



(33) 

mN, if a > 1 



C. On The Diversity Potential of Relay Selection 

Based on the high SNR analysis of the previous section, interesting results regarding the diversity order 
of relay selection can be obtained. These results are presented below, for different types of channels and 
different channel estimation techniques. 

1 ) Channel Estimation over Static Channels: Let us first focus on the scenario where the channel estimates 
are obtained via training in static channels. In this case, as can be seen from Table U the exponent a can 
take any positive value, depending on how fast the training power increases with SNR, with respect to the 
data transmission power. In particular, it follows from Theorem Q] that full diversity is achieved by using 
a training power which increases with SNR at least as fast as the data transmission power, i.e., a > 0. A 
slower increase with SNR results in a decreased diversity order. Further details regarding the latter argument 
are provided in Section IVIH via numerical examples. 
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2) FIR Channel Prediction in Time-Varying Channels: Here, we concentrate our attention on the case 
of FIR channel prediction; note that the special case of outdated channel estimates is also included in this 
scenario, by setting the number of predictor taps equal to one. As shown in Table HI this case results in a = 0. 
Interestingly, it follows from (l33T > that the diversity order equals the fading shape parameter, m, regardless 
of the number of available relays. In other words, when the estimates of time-varying channels are obtained 
via an FIR predictor, the diversity order of relay selection reduces to that of the scheme where only a single 
relay is available. From another viewpoint, when relay selection is performed over time-varying channels, 
its full diversity potential is completely lost, unless a predictor with an infinitely large length is employed. 
A study of the latter case follows. 

3) IIR Channel Prediction in Time-Varying Channels: As implied by Table U and Theorem [TJ the diversity 
loss of relay selection incurred by the time-varying nature of the underlying channels can be recovered via 
IIR channel prediction. Nevertheless, given that the full diversity order is recovered for a > 1, it follows that 
if the power of the channel estimation pilot symbols grows with SNR as fast as the data transmission power, 
i.e., a = 0, the resulting diversity order is still lower than the maximum value. As a result, it is concluded 
that in order to achieve full diversity in relay selection over time-varying channels, IIR channel prediction 
is required, in conjunction with a training power which increases faster than the data transmission power, 
i.e., a > 0. The amount of training power required to achieve full diversity is determined by the Doppler 
spread of the channel and the time difference between the consecutive noisy channel observations. 

VII. Numerical Results and Discussions 

Figs. Q] and [2] consider the case of static channels, and illustrate the effect of noisy channel estimates on 
the outage probability of relay selection over Nakagami-m fading. Specifically, Fig. Q] depicts results for the 
special case of Rayleigh fading, N = 3 available relays and (3 = 1, showing a significant dependence of 
the outage probability on the parameter a. Recall from Q that the parameters a and (3 reflect the relation 
between the power allocated to pilot symbols and the power used for data transmission. It is observed that 
full diversity is achieved for any a > 0, yet there exists a power gain loss compared to the perfect CSI 
case; this power gain loss is recovered for higher values of a, i.e., for a > 1. In Fig. \T\ we also observe 
the accuracy of the high SNR approximations in (l48l) . (l55l) . and (l50l) for a < 1 (a < 0), a « 1 (a « 0), and 
a > 1 (a > 0), respectively. 

The dependence of the outage probability on the ratio of the pilot power and data transmission power, 
(3, is depicted in Fig. [H for the case of static channels. Similarly as in Fig. [T] A = 3 available relays and 
Rayleigh fading (m = 1) are assumed, while the a parameter is set to a = 0. We notice a slight dependence 
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of the outage probability on f3 for a given a, which becomes negligible as /? grows large. Consequently, it 
is concluded that when the ratio of £ p and is constant in the whole SNR region, relay selection maintains 
its full diversity characteristics when operating over noisy static channels; the corresponding loss in power 
gain is noticeable only for small values of the ratio of £ p and 

Figs. [3] and |4] consider to the case of ideal but outdated CSI, where the channel estimation is assumed 
noise-free yet it suffers from feedback delay. Specifically, in Fig. [3] we assume the typical scenario of a 
vehicle moving at 50 km/h and receiving at a frequency of 2.4 GHz, which corresponds to a maximum 
Doppler frequency of approximately fd = 100 Hz. Under this assumption, we illustrate the dependence 
of the corresponding outage probability on the time interval between estimation updates, T d , for Rayleigh 
fading (m = 1), T = 1 and N = 5. 

We notice from Fig. [3] that the rate of estimation update significantly affects the outage performance of 
relay selection, in the sense that low update rates result in severe diversity and power gain losses. This is 
in agreement with (l33l where, given that for outdated channel estimates a = holds, the diversity order 
equals m regardless of T^. Nonetheless, it should be pointed out that for low values of Td the slope of 
the outage curves retains its full diversity characteristics in the practical SNR range, and approaches m 
only for infinitely high SNRs. This observation sheds light onto the diversity potential of relay selection 
with outdated channel estimates since, although it is impossible to achieve full diversity from a theoretical 
perspective (i.e., when j — > oo), it is still possible to achieve full diversity in the practical SNR range, by 
decreasing Fig. [3] also demonstrates that for relatively low channel estimation update rates (e.g., for 

= 3 msec), relay selection cannot take advantage of the large number of available relays, since the case 
of N = 5 yields approximately the same performance as that of no selection, i.e., N = 1. Furthermore, it 
is worth mentioning that the outage probability for the case where the mobile terminals are moving at the 
walking speed of 5 km/h can be also extracted from Fig. [3j by tenfolding the corresponding values of Td 
(i.e., T d = 30,25, msec). 

The outage probability dependence of relay selection with outdated CSI on the Nakagami-m parameter 
is depicted in Fig. 01 We assume five participating relays (N = 5) and an outage threshold SNR of T = 3, 
while the relation between fd and Td is set such that p = 0.5. As expected from (l33l . it is seen that increasing 
m results in a considerable outage probability decrease, accompanied by a shift of the slope of the outage 
curves at high SNR. Note that in all cases the diversity order equals m, as also corroborated by d33T ). 

On the basis of the moving-vehicle scenario considered above, which corresponds to fd = 100 Hz, the 
case of channel prediction in time-varying channels is treated in Figs. [5] and [6] In particular, Fig. |5]depicts the 
outage probability of relay selection in Rayleigh fading (m = 1) for several values of the channel predictor 
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length, L, including the case of L — > oo which corresponds to IIR channel prediction and serves here as 
benchmark. As demonstrated in Fig. [5j by increasing the number of predictor coefficients in FIR channel 
prediction the outage probability experiences a power gain increase, yet no diversity gain increase is seen for 
low values of L. On the contrary, an increase in the diversity gain is attained through IIR channel prediction, 
as shown in the ensuing, Fig. [6] 

Fig. [6] illustrates the outage probability of relay selection for IIR channel prediction, when operating over 
Rayleigh fading. We notice a high dependence of the outage probability on the value of Td, which corresponds 
to the time difference among the consecutive time instances in the infinite-length channel predictor. In 
particular, we notice that full diversity is achieved for small values of Td, while for larger values of Td the 
diversity characteristics of relay selection are lost, as expected from (|33T ) and Table U 

Finally, the achievable diversity order versus a for the unified imperfect CSI model, where the cases of 
noisy channel estimation and CSI imperfection due to time-varying channels are incorporated, is plotted in 
Fig. 13 As expected, we notice a linear increase of the diversity order for < a < 1 and constant diversity 
order for a > 1, which equals mN. 

In fact, Fig. |7] sheds an interesting light onto our general assessment regarding the diversity order of relay 
selection in Nakagami-m fading, which is as follows. The Nakagami-m fading model assumes multiple 
scatterers in each link, causing an "internal diversity" phenomenon of order m, which is independent of the 
channel estimation quality. On the other hand, the presence of multiple available relays offers the potential 
for additional, "external diversity" yet this additional diversity strongly depends on the quality of channel 
estimation, as reflected by a. Consequently, we notice from Fig. |3 that for any non-prime diversity order of 
mN, the value of m is more important than N for the overall diversity order for any or < 1; if a > 1, m 
and ./V affect the diversity order in exactly the same way. 

VIII. Conclusions 

We presented an assessment of the diversity potential of relay selection with practical channel estimation 
techniques, in Nakagami-m fading. The considered channel estimation techniques include the cases of 
estimation in noisy static channels; estimation in noiseless time-varying channels, and estimation in noisy 
time-varying channel with the aid of FIR and IIR channel prediction. A closed-form expression for the outage 
probability of relay selection with imperfect CSI was provided, as a function of the correlation coefficient, 
p, between the actual and the estimated channel values. Capitalizing on this outage expression, our principal 
inference was that the diversity order of relay selection is determined by the relative speed of convergence 
of p to one, compared to the speed that the SNR converges to infinity. 
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Appendix A 
IIR Channel Prediction 

Let us consider the process yi = hi + ni, where hi and ni denote the channel value and the corresponding 
noise component at time instance I, respectively, with I = 1, oo. The variance of the prediction error is 
derived as [25] 

/ rWT*) r / .„ ™x ?v„-| \ n„ 

(34) 



■l/(2Ti) 

a e = exp \T d I In 

-1/(2T„) 



5«i + ^ 



where Nq/£ p is the variance of the noise component and S^h (/) denotes the Fourier transform of (r). 
Therefore, considering that the spectrum of Shh (/) is band-limited by f d , d34l yields 



o-g = exp ( T d I In 



On 



d NoV~ 2fdTd _Ni 
Bp J £p 



o 



(35) 



Assuming fa > 0, (fTTT) is obtained from (I35T ). ©, and ©. 



Appendix B 

Derivation of the Auxiliary Function I (•,-,•,•, •) 

Using the alternative representation of the incomplete Gamma function shown in ll20l Eq. (8.352/4)], and 
applying the multinomial theorem, we obtain 

(TO— 1 ni i \ 3 



T J (m) exp (—j/3x) 



ni-l ((P x i 



71, 



(36) 



n o ,ni,...,n m _i=0 
no+ni+— +nm-i=j 



i=0 



Therefore, (1201 ) yields 

I(H,a,m,P,j) = 



E 

n o ,ni,...,n m _i=0 
n +ni+-.+n m _i=j 

P (m) r (j + 1) 

\M+1 



r» (m) r Cj + 1) I n r [f ! ^ 1 r / + S '\ 



i=0 



to-1 f £E 







7li ■ 



E in^r " + E 



(a + j/3f +1 ^ _ \ T(m + 1) 

y J r ' n o ,ni,...,n m _i=0 \ «=0 

n +rc,i+...+n m _i=j 

Appendix C 
Diversity Analysis in Nakagami-ttt, fading 



TO — 1 

// + > m,; I !. 

i=0 



(37) 



The following Lemma provides a high-SNR investigation of function G (m + j,a,b,l), allowing for a 
simplification of d3"TT ). 
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Lemma 3: For sufficiently high 7, the function G (m + j,a,b,l), defined in 02t . decays proportionally 

to fy-^olm+aj) t j_ e _ ? 

lim G(m + j, a,b,l) ~ . (38) 

7— >oo 

Proof: It follows from (1221) that the function ip (m + j,l — b^~ a , i) takes the following form 

if> (m + j, 1 - 67"°, i) = ^ + A % (m, 1 - 67"°, i) . (39) 

Therefore, for deriving the decay exponent of G (m + j, a, b, I) it suffices to evaluate the decay exponent of 
G (m,a,b,l). After substituting (1221) into 02t and algebraic manipulations, the function G (m,a,b,l) can 
be expressed as 

G (m, a, 6, = ^ — -— ^ j( } - .,. ^ (40) 

where £ G {0, | (Z — 1) (m — 1)}, are constants. It is observed from (l40b that the dominant term of 
G (m,a,b,l) in the high SNR regime decays in proportion to ^~ da , where cfe is given by 



c?G = am 



i J (to — 1) a + a [im — (i — 1)] = aZm. (41) 

\i=l ) i=2 



The proof then follows from (1411) and 091 , by applying the binomial expansion to the first term of the right 
hand side of 0JO). ■ 

For the derivation of the diversity order, we simplify ( f3TT ) for different values of a, as shown below. 

1 ) Case of a < 1: In this case, the second argument of T (to + n, mT/ (&7 1_a )) in (|3TT > tends to zero 
as 7 — ^ 00. Hence, simplifying the incomplete Gamma function similarly as in 001 ), we obtain 

/ mT \ m m+j y m+j 

T ( m+ ^- T { m + ^b^-a)^ (m+j)&m+ ^(l-a) (m+J r (42) 

Therefore, from 0TT) and (l42l ) we obtain an expression for the outage probability in the form of 

rwrO +C <43) 

where 

C = yV llll I V I V I G(rn+iabl) (44) 

j!r(m)r(m + j) (m + j)^^ 1 -*)^) Lr t m + * < WJ' ^ 

The following lemma investigates the high SNR behavior of C. 

Lemma 4: For sufficiently high 7 and a < 1, the quantity C defined in (1441 is approximated by reducing 

the sums to single terms corresponding to j = and / = iV, respectively, and decays in proportion to 

lim C ~ ^-m[a(AT-l)+l]_ (45) 
7— >oo 
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Proof: It follows from Lemma [3] that in the high SNR region the quantity inside the summation of (1441 
is analogous to j dc , where dc is given by 

dc = an — m (N — I) — (1 — a) (m + n) — (aim + an) = (1 — a) (ml — m — n) — mN. (46) 



Since a < 1, (1461 ) is maximized over the set of non-negative integers j for j = 0, regardless of I; this value 
of j corresponds thus to the dominant term in the outer sum of (l44l . Consequently, (l46l ) reduces to 

o!c = m [(1 - a) (I - 1) - iV] . (47) 

Since a < 1 and 1 < I < N, it follows from d4"71 ) that the dominant term in the inner sum of (1441) corresponds 



to Z = iV. Hence, setting I = N to (1471) completes the proof. ■ 
Therefore, considering the fact that the first term in (|43T ) decays in proportion to ^ _mAr ; it follows from 
Lemma |4] that for a < 1 and JV > 1, the dominant term in (|43T > is C, where only the term with j = and 
I = N is relevant, so that (14"3l reduces to 

rn m-lrpm 

P out « V(m) 6^7(1-^ G (m ' °' iV) ' (48) 
Eq. (l48l) represents a high SNR expression for the outage probability for a < 1. 

2) Case of a > 1: Let us now assume that a is larger than one, and does not lie in the proximity of one. 
The case where a approaches unity will be considered separately in Section |G31 For a > 1, we have 

mT \ 

T(m + j) . (49) 



Jim 

7— >oo 



r (m + j)-T[m + j 



Hence, following the same procedure as for proving Lemma 01 it is concluded that the Zth-order term within 
the double summation in d3TT ) decays proportionally to j~ m [ N + l { a - 1 )]^ irrespective of j. Consequently, since 
a > 1 the negative decay exponent of the second term of (1311 ) is higher than mN for any / > 1. This implies 
that the dominant term in (I3T1) is the first term for high SNRs, yielding 

/ m m-l\ N / T \mN 
Pout « - • (50) 



T (m) J \ 7 

Eq. (TSUI) represents the asymptotic outage expression in high SNR for a > 1. 

3 ) Case of a m I: The scenario where a lies in the proximity of one is treated separately, since in this 
case the second argument of T (m + j,mT/ (&7 1_a )) converges very slowly (to either zero or infinity) as 
7 — > 00, hence the approximations in (l42l ) and (|49l do not hold for practical SNR values. As a result, the 
transition from the a < 1 case to the a > 1 case experiences a discontinuity in terms of the (practical) high 
SNR approximation of the outage probability, as a approaches unity. This discontinuity is bridged through 
the outage expression presented below. Recall from Section IIII-BI that the case of a = 1 corresponds to the 
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common scenario of channel estimation in noisy static channels, where the power allocated to pilot symbols 
equals the power allocated to data transmission. 

Since awl, let us assume that 7 1 ~ a approaches a non-zero finite constant as 7 — > 00, i.e., 

lim f~ a = A, < A < 00. (51) 

7— >oo 

This allows us to evaluate the integral I (m + j — l,m/ (&7 1-a ) ,m,m/j,i), shown in (l20l) . as follows. 
• Let us assume < x < 00. Then, using [ 20 , Eq. (8.352/7)] we have for high SNR 



r (m, yxj « P (m) . (52) 

Let x — >• 00. In this case, (l52l) does not hold. However, it follows from LHospital's rule, as well as 
from the fact that bX is finite, that 



lim 

x— >oo 



x^^expi—x) V m,-x 

OA J \ 7 



= 0. (53) 

V 7 /J 

Therefore, it follows from (1521) and ([53]) that for high SNRs 

I (m + j - l,^^,m,™,i \ ^ x m+j - 1 ^ V ( K -^jY i {m)dx = T i {m)T{m+j) {^) 

(54) 

By combining (fJD, (EH), 63, <EB, <E3>, lEOl Eq. (0.15.4)], [20, Eq. (8.310/1)], E3 Eq. (8.350/2)], and 



the infinite series representation of the exponential function, we arrive after some manipulations at 



,m-l\ N / m \ mN 



Pout^(N + l)[^—) (?) . (55) 



T ym) J \ 7 

Eq. (|55T ) represents a high SNR approximation of the outage probability for the case where a lies in the 
neighborhood of one. Theorem Q] follows then directly from (|48l , d55l ), and (|50l ). 
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TABLE I 

Parameters a and b for the considered imperfect CSI scenarios. 



Case 


a 


b 


Outdated CSI 





1-P 


Noisy CSI 


a + 1 


ml {PL) 


FIR Channel Prediction 







IIR Channel Prediction 


(a + 1) (1 - 2f d T d ) 


exp (t„ /f- d In (e^ 2 ^^)] df) p-W*™ 




5 10 15 20 25 30 35 40 

Average SNR y (dB) 



Fig. 1. Case of static channels: Outage probability of relay selection with noise-impaired channel estimates in Rayleigh fading 
(m = 1) versus average SNR, assuming N = 3 available relays, T = 1, /? = 1 and several different values of a. 
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Fig. 2. Case of static channels: Outage probability of relay selection with noise-impaired channel estimates in Rayleigh fading 
versus average SNR, assuming N = 3 available relays, T = 1, a = and several different values of j3. 
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Fig. 3. Case of outdated CSI: Outage probability of relay selection in Rayleigh fading versus average SNR, assuming N = 5 
relays, T = 1, maximum Doppler frequency fa = 100 Hz and several different values of Td. 
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Fig. 4. Case of outdated CSI: Outage probability versus average SNR for p — 0.5, N = 5, T = 3, and several values of the 
Nakagami-m shape distribution parameter, m. 



August 22, 201 1 



DRAFT 



25 




Fig. 5. Time-varying channels with FIR channel prediction: Outage probability of relay selection in Rayleigh (m = 1) fading 
versus average SNR, assuming N = 5 available relays, T = I, T d — 3 msec, f d = 100 Hz, a = 0, /3 = 1, and several values of 
the channel predictor length, L. 



August 22, 2011 



DRAFT 



26 



CO 
_Q 
O 



0.01 
1E-3 
1E-4 
1E-5 
1E-6 
1E-7 
1E-8 
1E-9 



* -4- v . • 



CD 

CD 

CO 
-*—> 

° 1E-10 



1E-11 
1E-12 
1E-13 



1 □ 
l-O- 

-A- 

^--v- 



T =5msec\ s n% 
T =4msec 

d 

T =3msec 

a 

T =2msec 

a 

T = 1msec 

d 

T =0. 5msec 

d 



;— -0--T. =0.1 msec 



5 



~ l 

10 



1 5 



I 

20 



25 



~~ r 

30 



V' 



"a. 



35 



40 



Average SNR y (dB) 



Fig. 6. Time-varying channels with IIR channel prediction: Outage probability of relay selection in Rayleigh (m = 1) fading 
versus average SNR, assuming N = 5 available relays, T = 1, a = 0, /3 = 1, and several values of Ta- 



August 22, 2011 



DRAFT 



27 




Fig. 7. Diversity order of relay selection with imperfect CSI in Nakagami-m fading as a function of the exponent parameter a, 
for several combinations of the number of available relays, N, and the fading shape parameter, m. Applies to all the considered 
imperfect CSI scenarios. 
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